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Abstract. Gamma-ray binaries are suitable sources to study high-energy processes in jets 
and outflows in general. In the last years, there has been a lot of activity in the field of 
gamma-ray binaries to identify the different factors that shape their non-thermal spectra, 
which ranges from radio to very high energies, as well as their lightcurves. In this work, I 
discuss the main aspects of the non-thermal emission in this class of objects, which presently 
includes high-mass microquasars, high-mass binaries hosting a non-accreting pulsar and, 
probably, massive star binaries; few potential candidates to be gamma-ray binaries are also 
presented. Finally, the importance of gamma-ray absorption is discussed, and the main phys- 
ical ingredients, which are likely involved in the non- thermal radiation in gamma-ray bina- 
ries, are briefly considered. 
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1. Introduction 

Gamma-ray binaries are binary systems 
formed by two objects, generally one non- 
degenerated star and a compact object or 
another non-degenerated star, which emit 
gamma-rays. 

Three gamma-ray binaries have been 
detected so far in high- (HE: ~ 0.1 - 100 GeV) 
and very high-energy (VHE: ~ 0.1 - 100 TeV) 
gamma-rays: LS 5039, LS I +61 303 and 
Cygnus X-l ( at ~ 4 cr) dKniffen et al l 1 19971: 
Pare des et all 2000; Aharonia n'etaD |2()05a, 
2006 d: lAlbert et al l 12006. 2007: lAcciari et alJ 
2008; lAbdoetalJ l2009allbl: ISabatini et al.l 
2010). In addition to these three sources, there 



Send offprint requests to: V. Bosch-Ramon 



are a gamma-ray binary and a candidate for 
it, PSR B 1259-63 and HESS J0632+057, 
respect ively, detected above 100 GeV , but not 
below (lAharonian et alJl2005bll2007b . Finally, 
Cygnus X-3 has been de tected in the HE but 
not in the VHE band dTavani et all l2009al: 
lAbdo et al.1 l2009d: ISaito et al.l 120091) . as well 
as probably 77-Carina, which cou ld have been 
detecte d only in the HE range dTavani et al.l 
2009b). Among all these sources, there are two 
confirmed high-mass microquasars (HMMQ) 
(Cygnus X-l and Cygnus X-3), one confirmed 
high-mass binary hosting a non-accreting 
pulsar (PSR B 1259-63), possibly a massive 
star binary (77-Carina), two high-mass binaries 
harboring a compact object of unknown n ature 
(LS 5039 and LS I +61 303; see, e.g. JPubus! 
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2006at IChernvakova et al]l2006URomero et al.1 
120071: iBosch-Ramon & Khangulyanl 

2009), and a high-mass binary candidate 
(HESS J0632+057). 

In the following, we briefly introduce the 
different types of source that have been de- 
tected in gamma-rays (Sect. [2]), as well as oth- 
ers that may be detectable in the nearby fu- 
ture (Sect. [3]). The importance of gamma-ray 
absorption in gamma-ray binaries is also dis- 
cussed (Sect.©, and a short inventory is made, 
together with some remarks, of all the differ- 
ent physical ingredients that are probably to 
be considered in these sources (Sect. [5]) . For a 
very recent observational overview of gamma- 
ray b inaries, we address the reader to IParedei 
(20l3). 



2. Physical scenario of established 
gamma-ray sources 

2.1. High-mass microquasars 

Microquas ars are X-ray binaries with rad i o jets 
(see e.g., iMirabel & Rodriguezl Il999h iRibol 
I2005I) . and when hosting a massive star, they 
are called HMMQ. The power supply in these 
objects can be either accretion or an accret- 
ing rotating black-hole, and emission would 
be directly powered by a jet launched in the 
inner regions of the accretion disk. The jet 
can produce non-thermal populations of rela- 
tivistic particles via diffusive shock accelera- 
tion or other mechanisms at di fferent spatial 
scales (for a recent rev iew, see [ Bosch-Ramon 
201 Oat see also, e.g., iLeyinson & Blandford 



1996 :lRomero et al.ll2:QQ3l: iParedes et al. 2006; 
Rieger et al.ll2007l) . In the context of HMMQ, 
the most efficient gamma-ray mechanism 
would be inverse Compton (I C) scattering 
of the stellar photons 



Paredes et al 



2000t iDermer et al.ll2QQ6L IBosch -Ramon et al 
2006al) . A sketch of the microquasar scenario 
is shown in Fig.[T] 

Beside the jet itself, stellar wind-jet dy- 
namical interactions at the binary spatial 
scales could also lead to non-thermal emis- 
sion in HMMQs (iPerucho et al.l 120081) . and 
that the likely clumpy nature of the wind 



(lOwocki & Co hen 2006) could lead to HE and 
VHE flares (see, e.g jAraudo et al.ll2009l) . 

An example of the importance of the wind- 
jet interaction is shown in Fig. [2 in which it 
is shown a 2-dimensional map of the density 
of a jet of power 10 35 erg s" 1 that is launched 
and interacts with the stellar wind, which is 
coming from the top. Note that the jet is suf- 
fering significant bending and strong disrup- 
tion not far beyond an asymmetric recollima- 
tion shock. Jets significantly more powerful, of 
up to 10 37 erg s" 1 , may be still affected by jet 
disruption due to the growth of non-linear in- 
stabil ities triggered in t he recollimation shock 
fsee lPerucho et al.ll20ldh . 

2.2. High-mass binaries with a 
non-accreting pulsar 

In young pulsar binary systems, the non- 
thermal emission is expected to be generated 
in the region of collision between the mas- 
sive star (non-relativistic) and the pulsar 
(relat ivi stic) winds (e.g . [Maraschi & Treves 
19811: iTavani & Aronsl Il997l: iKirk et al 



119991: iDubusI 120063 [Khangulvan et al 



120071: iNeronov & Chernvakoval |2007), al- 
though the shocked pulsar wind, advected 
away and accelerated by s trong pressure 
gradients (iBogovalov et al.l [2008), could 
also b e a source of Do p pler bo osted emis- 
sion dKhangulyan et al.l l2008al) . The un- 
shocked pulsar win d may be also a source o f 
gamma-rays (e.g. iKhangulyan et al.l [2007), 
which may suffer reprocessing in the stel- 
lar field due to pair creat io n (se e, e.g., 
Sierpowska-B artosik & Torre"sl 120071) . This 
reaccelerated flow may release kinetic energy 
far from the shock region, increasing the 
complexity of the phenomenology of the 
standard scenario, since this could affect the 
non-thermal emission in the whole spectrum. 
A picture representing the scenario discussed 
in this section is shown in Fig. [3] 



2.3. Massive star binaries 

Massive star binaries are systems in which a 
strong shock takes place between the winds 
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Fig. 1. Illustrative picture of the microquasar scenario, in which the main elements and processes relevant 
to this work are shown (background image adapted from ESA, NASA, and Felix Mirabel). 



Logarithm of rest-mass density. X cut 




Logarithm of rest-mass density. Z cut 




Fig. 2. Two-dimensional map of the density of a jet interacting with the stellar wind, which is coming from 
the top. 



of the stars. These shocks have speeds of 
up to few 1000 km s" 1 and are collision- 
less and strongly supersonic, and particle ac- 
celeration and non-thermal emission up to 
gamma-ray energies, m ainly through IC , could 



be efficient in ther e (lEichler & Usovl 



1993 



Benaglia & Romerol [20031; iBednarekT I2OO5I: 
Rei mer et alJ 120061: lDe~Beckerl 120071: iPittardl 
20101) . In Fig. g] we show the result of a 2- 



dimensional simulation of the hydrodynamical 
interaction between the winds of two O stars. 
The strong shock between the two stellar winds 
is clearly visible in the density jump. 



3. Other potential gamma-ray binaries 

Beside HMMQ, non-accreting pulsar high- 
mass binaries and massive star binary systems, 
all of them already detected in gamma-rays, 
other binaries are potential candidates to be 
also found in this energy range. 

Low-mass microquasars have been also 
proposed to be HE and VHE emitters. In 
these objects, the most efficient gamma-ray 
mechanisms could be external IC with accre- 
tion photons, synchrotron self-Compton or 
hadronic interactio ns in the inner-most regions 
of the jet base (e.g.|Atovan & Ah aronian 1999; 
Bosch-Ramon et al.1 l2006at iRomero & Vilal 
20091) . [n microquasars in general, the corona 
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Fig. 4. Two dimensional map of the density distribution in a system of two O stars with colliding winds. 
The wind of the star to the left is 10 times denser than that of the star to the right. 




Fig. 3. Illustration of a binary system consisting of 
a pulsar and a Be star. The pulsar orbits a massive 
Be star with a disk- like outflow of stellar material 
(from Hubble archive). Around periastron, the pul- 
sar outflow interaction with the stellar wind is the 
strongest, leading to particle acceleration and non- 
thermal emission. More regular and circular systems 
may yield stable emission all along the orbit. 



(or the base of the jet) could be also a non- 
therm al emitting region (e.g. iGierlinski et al.l 
119991: see also IRomero et al.1 l20ld) . as well 
as th e termination reg ion of the jet in the 
ISM (iBordas et al.ll2QQ9h . In Fig. S we show 
the computed spectral energy distribution for 
a powerful jet of a low-mass microquasar. 
The dominant IC component is either IC 
scattering off corona photons or synchrotron 
self-Compton; note the soft spectrum above 



rona/sync 

jet dorr 
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ORXTE 
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Fig. 5. Computed spectral energy distribution of the 
non-thermal emission from IE 1740.7-2942 for two 
situations. In one case, the hard X-rays come from 
a corona (long dashed line), whereas in the other, 
they are of synchrotron origin and come from the jet. 
Absorption in the accretion disk and corona photon 
fields has been taken into account (see the two dips 
below and above ~ GeV energies). For details, see 
Bosch-Ramon et al. (2006b). 



GeV energies due to the Klein Nishina effect. 
Gamma-ray absorption in the accretion disk 
and corona fields can b e significant in this 
contex t (for details, see iBosch-Ramon et all 
I2006bl) . 

Accreting pulsar X-ray binaries might be 
also sources of ga mma-rays, as i t has been 
suggested by, e.g., IRomero et al. I (1200 ll) and 
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Sgue ra et al.l (120091) . each of these works point- 
ing to different powering mechanisms: pro- 
ton beams moving along the pulsar magneto- 
spheric magnetic lines and colliding with the 
accretion disk, or sudden ejections of matter 
through the formation of a magnetic tower pre- 
ceded by high-rate accretion events. This high- 
rate accretion events would lead to overcome 
the magnetic pressure of the pulsar magneto- 
sphere allowing for jet formation (see, e.g., 
M assi & Kau fman Bernado 2008, and refer- 
ences therein). 

4. Impact of gamma-ray absorption 

Other non-thermal emitting regions dif- 
ferent from those already discussed 
could exist in any type of compact high- 
mass binary. If gamma-ray absorption 
takes plac e under th e photon field of 



the star dFordl 1 19841: 



Protheore & Stanevl 
19871: iMoskalenko & Karakul j 
Boettcher & Permed 120051: iDubusl 120061 
Kha ngulvan et al.l l2008bt iRevnoso et all 
2008b . the whole binary system could be 
an efficient broadband non-thermal emitter. 
In such a case, the created pairs can radiate 
through synchrotron and IC emission in the 
whole spectral band, interacting with the 
ambient mag netic field and stellar ph otons, 
respectively (iBosch-Ramon et ap2008al) . 

In case the energy density of the mag- 
netic field is much smaller than that of the 
stellar photon field, IC becomes the dominant 
cooling channel of these pairs. Then these 
pairs can produce more gamma-rays that to 
their turn may be also absorbed, triggering 
what is called an electromagnetic cascade. 
In this way, the effective optical depth to 
gamma-rays can be significantly reduced. The 
deflection of the created pairs in the ambient 
magnetic field determines whether t he cascade 
is lin e ar or three dime nsional (e.g. Bed narekl 
1997; Aharo nian et al. 



2007 



... I2006bt lOrellana et al I 
Sierpowska-Bartosik & Torresl 12007 ; 



Zdziarski et al.ll2009l:ICerutti et alJl2010h . 

If the magnetic field is high enough, 
electromagnetic c ascades are suppressed 
dKhangulvan et al.l l2008bh and the X-ray 
(synchrotron) and GeV (single-scattering 
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Fig. 6. Computed spectral energy distribution of the 
radiation produced by secondary pairs in LS 5039 
(green lines) at the orbital phase associated to the 
superior conjunction of the compact object (the bi- 
nary system properties can be found in Aragona et 
al. 2009). Two different configurations have been 
taken, one in which the emitter is at a height of 
10 12 cm (production gamma-ray spectrum: thick 
solid blue line; corresponding pair injection: thick 
solid red line), and at 2 X 10 12 cm (production 
gamma-ray spectrum: thin solid blue line; corre- 
sponding pair injection: thin solid red line) above 
the compact object. The production spectral energy 
distribution (thin long-dashed) of the secondary IC 
emission is also shown. 



IC) emission from secondary pairs can 
dominate the secondary energy output 
(Bosch-Ramo n et al . 2008a). As two examples 
of this, we show in Figs. [6] and [7] the spectral 
energy distribution of the secondary radiation 
for a gamma-ray binary with system properties 
similar to those of LS 5039. The magnetic 
field in the stellar surface has been fixed to 
20 G. For stronger enough magnetic fields, 
the secondar y radio emission may be also 
significant dBosch-Ramonl I2009L 1201 Obh . In 
Fig. [8] we show 5 GHz radio images for a 
similar case but adopting a star magnetic 
field of 200 G. The images for four different 
orbital phases are presented: 0.0 (superior 
conjunction of the compact object), 0.25, 0.5 
(inferior conjunction), and 0.75. 
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Fig. 8. Computed image of the 5 GHz radio emission, in the direction to the observer, for different 
orbital phases. Units are given in mJy per beam, being the beam size ~ 1 milliarc second. 
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Fig. 7. The same as in Fig. [6] but for the inferior 
conjunction of the compact object. 



5. Main elements to understand the 
gamma-ray emission 



The phenomenology at high energies of 
gamma-ray binaries can be (at least) semi- 
quantitatively explained accounting for a set 
of ingredients or physical processes and ef- 
fects: the angular dependence of gamma-ray 
absorption and IC scatteri ng and the (mod- 
erate) role of absorption ( Khangu lyan et al.l 
2008b), in some cases locating the emit- 
ter at some height above the orbital plane 
(for acceleration and mode ling arguments, 
see iKhangulyan et al.l l2008bt for absorptio n 
arguments, see iBosch-Ramon et al.l [2008b), 
likely accounting for adiabatic losses (e.g . 
Khangulvan et al. 2007; Takahashi et al. 2009; 
Zaba lza et al.ll2010b . The role of the magnetic 
field, e.g. through the impact of synchrotron 
cooling, is crucial for the VHE spectrum and 
luminosity. In HMMQs, one should not ne- 
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gleet the impact of the stellar wind on t he non- 
thermal processes dPerucho et al .1120 101) . 

All the mentioned physical effects and pro- 
cesses are present to different extents in all 
the known types of gamma-ray binaries. This 
makes in some cases the identification of the 
specific nature of the object a difficult task, 
as it is shown by the debate on the nature of 
LS 5039 and LS I +61 303, i.e. HMMQs ver- 
sus non-accreting pulsar high-mass binaries, if 
not a really new class of object. Therefore, it 
is of primary importance to apply comprehen- 
sive approaches to study theoretically these ob- 
jects, including semi-analytical modeling, nu- 
merical calculations (secondary emission) and 
(magneto)hydrodynamical simulations, as well 
as carry out simultaneous multiwavelength ob- 
servations with high angular, spectral and time 
resolution. Nowadays, the field of gamma-ray 
binaries seems to be the best example for the 
need of such a complex multidisciplinary ap- 
proach, if one wants to understand the details 
of the physical processes taking place in these 
objects, as well as to distinguish between dif- 
ferent scenarios and models. 
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